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Abstract. Helical multishell gold nanowire is studied theoretically for the formation 
mechanism of helical domain boundary. Nanowires with the wire length of more 
than 10 nm are relaxed by quantum mechanical molecular dynamics simulation with 
tight-binding form Hamiltonian. In results, non-helical nanowires are transformed into 
helical ones with the formation of atom pair defects at domain boundary, where the 
defective atom pair is moved from an inner shell. Analysis of local electronic structure 
shows a competitive feature of the energy gain of reconstruction on wire surface and 
the energy loss of the defect formation. A simple energy scaling theory gives a general 
discussion of domain boundary formation. 
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1. Introduction 

Nanometer scale material forms exotic structures and it is crucial to understand and 
control their formation mechanism, so as to establish the foundation of nano electronics. 
Helical multishell gold nanowire [H [2] is one of exotic metal nanostructures [3j and has 
characteristic multishell configurations with 'magic numbers'. The helical wires were 
synthesized by focusing electron beam on thin film. [1] The wire axis is in the [110] 
direction of the original FCC structure and the outermost shell is a folded (lll)-type 
(hexagonal) atomic sheet with helicity. A single shell helical structure was synthesized 
later. [5] The transport property was studied theoretically [6] and experimentally [Tj. 
Platinum nanowires were also synthesized with the same type of helicity. [8] 

The observed multishell structures of the helical gold nanowires [1] are denoted 
by the numbers of atoms in each shell and are classified into 7-1, 11-4, 13-6, 14-7-1, 
15-8-1 structures. For example, a '14-7-1 nanowire' is a nanowire with three shells and 
the outer, middle and inner shells have fourteen, seven and one atom(s) in the section 
view, respectively. These numbers are called 'magic numbers', since the difference of 
the numbers between the outermost and next outermost shells is seven, except the cases 
of the 7-1 structure. 

Theoretical investigations were carried out for the helical nanowire structures 
P[IOl[ni[T2llHll3l[l5l[T6llI7] and, among them, a recent theory ^ [IE] explains 
the observed multishell configuration with 'magic numbers' systematically. The theory 
was proposed as a two-stage formation model of helical multishell gold nanowires and 
the model was confirmed by quantum mechanical molecular dynamics simulations for 
gold and copper with tight-binding form Hamiltonians. The simulation shows the 
energy relaxation process from non-helical structures into helical ones. The analysis 
of electronic structure concluded that the origin of the helicity comes from the intrinsic 
nature of non-spherical 5d electrons and a (lll)-type (hexagonal) surface structure is 
energetically favorable for 5d electrons in sheet structure. The above mechanism gives a 
general understanding among (i) the appearance of helical nanowire structures of gold 
and platinum and (ii) the fact that reconstructed equilibrium surfaces of FCC 5d metals, 
gold platinum and iridium, prefer (lll)-type structures. [HI [201 [2ll [221 [231 [211 [25] After 
the theory paper, [15] several related simulations were carried out for formation of helical 
gold nanowire within tight-binding form Hamiltonian. [161 HI] 

The present paper investigates the formation of helical domain boundary on wire 
surface, where the fundamental picture and simulation method are shared with the 
previous theory paper. p31| In general, a defect should be introduced at domain 
boundary in the formation process of a helical domain from non-helical one. The 
information of domain boundary is missing in the previous paper, since the simulated 
structures are short isolated wires of which wire length is less than 3 nm. In the 
resultant helical nanowires, the domain boundaries are located in the wire ends that 
are terminated artificially. In the present paper, simulations are carried out with longer 
wires, longer than 10 nm, and reveal a possible defect induced mechanism for forming 
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Figure 1. Simulated formation process of helical multishell gold nanowires with 
the wire lengths longer than 10 nm. The 11-4 structure is depicted in (a) at the 
initial structure (left panel) and the 9,000-th iteration step (right panel) and the 15-8- 
1 structure is in (b) at the initial structure (left panel) and the 6,600-th iteration step 
(right panel). See text for details. 

helical domain boundary. 

This paper is organized as follows; Section [2] describes the methodology and result 
of the simulation. Section [3] focuses on the analysis of the results, particularly on the 
energy per atom and the local density of states (LDOS). In Sec. HI a simple energy 
scaling theory with respect to the wire length of helical domain is constructed for a 
general mechanism of domain boundary formation. Finally, the summary and future 
aspects are given in Sec. [H 

2. Simulation 

The formation process of the helical wires is realized by quantum mechanical 
molecular dynamics simulation as a structure relaxation with thermal fluctuation. The 
temperature was set to be T =600 K, lower than the melting temperature (1337K). 
The simulation code used here has the name of 'ELSES' (=Extra-Large-Scale Electronic 
Structure calculations) . |26l|27l[28l[29l|30l[ni[33[33l[T5l[I^ simulation was 

reahzed with a tight-binding form Hamiltonian [36l [371 EH], which was used for gold 
nanowire. [iQl HH [391 [IS IIHl [13 The Hamiltonian form contains several parameters 
and they are determined to represent electronic structures of bulk solids, surfaces, 
stacking faults and point defects. One iteration step corresponds to At = 1 fs. Eigen 
state calculation for electronic structure is carried out at each iteration step. The 
boundary condition is imposed by fixing the mass center of the wire end layers. [T5] 
The computational time with 608 atoms is typically nine minutes per one iteration step, 
when a workstation with four dual-core Opteron^"""^) processors was used. 
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We should say that the present simulation is different from experiment in several 
points. For example, the simulated time is quite short, 7-9 ps, owing to the practical 
limit of computational resource. Therefore, the simulation result should be understood 
so that it captures an intrinsic energetical mechanism of the real process. 

The resultant nanowires are shown in Fig.[T]for the (a) 11-4 and (b) 15-8-1 structures 
and contain multiple helical domains on wire surface with well-defined domain boundary, 
unlike the smaller samples in the previous paper. [15] The initial structures are parts 
of ideal FCC structure and are the same as those in the previous paper, except their 
wire length. The wire lengh of these nanowires is approximately 12 nm or is 40 or 
42 unit layers for the 11-4 and 15-8-1 structures, respectively. Here a unit layer of a 
[110] nanowire consists of two successive atom layers. At the initial structures, the wire 
surface consists of (lll)-type (hexagonal) and (OOl)-type (square) areas. The latter area 
consists of two atom rows that are indicated by the two parallel arrows in Fig. [H Defect 
structures appear on wire surface as helical domain boundary. Hereafter the domain 
boundary formation of the 11-4 structure will be focused on. 

Figure [2] shows several snapshots of sideview for the 11-4 case, where the atoms are 
distinguished by color, so as to clarify the transformation process. Red atoms are those 
initially placed on the (OOl)-type surface area and green atoms are those initially placed 
in the inner shell region. The other atoms are those initially placed on the (lll)-type 
surface area. The central region indicated by broken lines at Fig. [2]|^a)-(d) is sketched 
in Fig. [2](e) (f) and (g). Figure [2t^e) and (f) capture the structures of Fig. M^a) and 
Fig. [2](b), respectively and Fig. [2] (g) captures the structure of Fig. [2]^c) or Fig. EJ^d). 
At the initial structure. Fig. [2](e), the 'red' atoms, marked as P, Q, V, W, X and Y, form 
a square lattice on wire surface and the 'green' atoms, marked as R and S, are placed 
in the inner shell. At the final snapshot. Fig. [2](g), the two 'green' atoms R and S, are 
moved from the inner shell into the wire surface, which form a helical domain boundary. 
They are inserted between the atoms P and V, while the atom Q is moved into an upper 
area between the atoms T and U. As result, the surface reconstruction occurs, from 
the (OOl)-tj^e (square) lattice into the (lll)-type (hexagonal) one, and introduces the 
helical domains with domain boundaries. 

As a remarkable tendency of the entire sample, surface defects at domain boundaries 
appear typically as pairs of 'green' atoms, as illustrated in Fig. MJ^) and (i). The 
central atom-pair defect, the pair of the atom R and S, appears, because the shear-like 
deformation on the (OOl)-type lattice occurs in the opposite shear directions between 
the left and right areas of the defect, as indicated by a blue arrow in Fig. EJ^i). 

3. Analysis of local electronic structure 

The transformation process from the non-helical wire into the helical one occurs with 
energy gain, according to the change of the total energy shown in Fig. [3](a). 

The change of the atom energy, energy of each atom, is plotted in Fig. [3]^b) for the 
selected atoms P, Q, R and S, that are indicated in Figs. [2] (e), (f) and (g). Figure IHa) 
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Figure 2. The transformation process of the 11-4 structure, (a) (b) (c) (d) : 
Snapshots at the initial structure and 1,000-th, 3,000-th and 9,000-th iteration steps, 
respectively, (e) (f ) (g) : Sketched figure as a close up of the central area that are 
indicated by broken lines on the wire surface. The figures (e) and (f) capture the 
structures of (a) and (b), respectively and the figure (g) captures the structure of (c) 
or (d). The atoms depicted as smaller balls, R and S in (e) and S in (f), are placed 
within the inner shell, (h) (i) : Schematic figure of formation of atom-pair defects in 
the (h) initial and (i) final structures. 

shows the LDOS of the atom P at the initial structure and the 1,000-th and 9,000-th 
steps. Figure Hl^b) shows the LDOS of the atom R at the initial structure and the 
1,000-th, 1,300-th and 9,000-th steps. In Fig. HI the origin of the vertical axis is shifted 
among the snapshots and the origin of the horizontal axis is set to the Fermi level for 
each snapshot. Since the values of the Fermi level are different only about or within 
0.1 eV among the snapshots, the difference in the Fermi level is negligible in the energy 
scale of the graphs. Since the atoms R is placed in the inner shell at the initial structure, 
the LDOS profile, shown in the lowest graph of Fig. 111(a), is similar to that of bulk FCC 
gold and has characteristic three peaks in the 5d band at £ ~ — 2 eV, -4 eV and -6 eV. 
The atom movement from the inner shell into the wire surface accompanies the 
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Figure 3. Energy change of the 11-4 structure during the relaxation process, (a) 
Change of the total energy, (b) Change of the atom energy for the atom P, Q, R and 
S that arc indicated in Fig. [2l 
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Figure 4. Local density of states of the atoms P and R at the initial structure and 
several iteration steps written in the graphs. The origin of the vertical axis is shifted 
among different iteration steps. The origin of the horizontal axis is set to the Fermi 
level. 
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drastic change of the atom energy. Figure [21(b) contains temporal peaks nearly at the 
1,000-th and 2,400-th iterations for the atom R and S, respectively, which corresponds 
to the successive atom movements of the atom R and S, from the inner shell into the 
wire surface, shown in Fig. [2l 

When the initial and final structures are compared in the atom energy, one can find 
that the helical transformation is motivated by a surface effect. At the final iteration 
step (9,000-th iteration step), the four atoms P, Q, R and S are transformed into 
members of the folded (lll)-type surface with helicity, as shown in Figs. [2] (d) and (g). 
Therefore the atom energy of them reaches an almost unique value, E ^ —3.1 eV in 
Fig. Mjo), and the final LDOS profile of the atoms P and R are quite similar, as plotted 
in the highest graphs of Fig. |l](a) and (b). In Fig. [2l the atoms P and Q, depicted as 
'red' atoms are transformed from the (OOl)-type surface area into the (lll)-type one. 
Figure [3](b) shows that the atom energy of the two atoms decreases from E ^ —2.5 eV 
into E f=^ —3.1 eV. The energy gain mechanism is explained by the two-stage model in 
the previous paper and the present case corresponds to the case of the atom B in the 
paper. [15] On the other hand, the atoms R and S, depicted as 'green' atoms in Fig. [2l 
are moved from the inner shell into the (lll)-type surface area. Figure [3]^b) shows that 
the atom energy of the two atoms increases from E ^ —4.3 eV into E ^ —3.1 eV. 

The above observation in the increase and decrease of the atom energy is confirmed, 
when histograms (not shown) are constructed for the atom energy difference AE between 
the initial and final structures. The histograms are constructed among all the atoms 
except the 158 atoms located within the five layers near the two wire ends, so as to avoid 
possible artifacts of the wire ends. All the 'red' atoms, the atoms originally placed on the 
(OOl)-type area of the surface, show the energy decrease of AE ^ —0.2 ~ —1.2 eV and 
the histogram peak is located at AE ^ —0.7 eV. The 'green' atoms, the atoms originally 
placed in the inner shell, show the energy difference from AE ^ - 0.2 eV to 1.6 eV. The 
large energy increase {AE 1.2 ~ 1.6 eV) comes from the atoms moved into the wire 
surface, like the atom R and S. Energy increase also appear among several other 'green' 
atoms in the inner shell, since the resultant helical wire has defective regions in the inner 
shell, particularly regions near the moved atoms. The histogram of the rest atoms, the 
atoms originally placed on the (lll)-type area of the surface, shows a wide range of 
values AE —1.2 ~ 1.0 eV and their average is an energy decrease {AE ^ —0.16 eV). 
We should remember that these 'rest' atoms are different in their situations in the initial 
structure. See the cases of the atom A and C in the previous paper. [15] 

It should be noted that, an experimental paper of [110] gold nanowire earlier 
than the report of helical structurefl], suggests a surface stabilization mechanism in 
which (OOl)-type area on the wire surface reconstruct into (lll)-type ones as on 
equilibrium surface. The suggested stabilization mechanism is consistent to the previous 
theory [15j and the present analysis. 

The analysis of the intermediate structures reveals a competitive feature of temporal 
energy gain and loss among atoms. In Fig. [3](b), for example, the atom R shows energy 
loss during the 800-th and 1,000-th iteration step {E f» —3.7 eV =^ —2.7 eV), while 
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the atoms Pand Q show energy gain during the same period {E ^ —3.1 eV =^ —3.7 
eV and E ^ —2.8 eV =^ —3.4 eV for P and Q respectively). The temporal energy gain 
and loss are almost the same amplitude. During the period, the atom R is moved from 
the inner shell to the wire surface (See Fig. [2] (b)). In conclusion, the atom R is moved 
with energy loss, from the inner shell into the wire surface, because of the energy gain 
mechanism of the atom P and R. The energy gain is found also in the LDOS profile 
of the atom P, shown in Fig. IH^a). When the LDOS profile is compared between those 
at the initial structure and at the 1,000-th step, the weighted center of the 5d band is 
shifted downwards. Moreover a gap-like structure appears, at e ~ —4.5 eV, in the LDOS 
profile of the 1,000-th step, which indicates the formation of bonding and antibonding 
states. In short, an unstable surface atom is changed into a ( relatively ) stable one, 
owing to the energy gain of the 5d electrons. [I5] The atom R also shows a temporal 
energy gain after the movement into the wire surface, during the period between 1,000- 
th and 1,300-th iterations (E f» -2.7 eV -3.8 eV). The change of the LDOS profile 
during the period, shown in Fig. Hl^b) is similar to the change of the atom P discussed 
above, which indicates the same stabilization mechanism with the 5d electrons. 

4. Discussion 

The above analysis indicates that the helical transformation is motivated by the energy 
gain in the surface reconstruction from the (OOl)-type area the into (lll)-type one and 
sacrifices the energy loss for forming point defects in the domain boundaries on the wire 
surface. 

A simple theory is constructed by energy scaling with respect to the wire length 
of a helical domain Ldom- The main energy gain appears among 'red' atoms that are 
distributed on line and the energy gain is scaled as O(Ldom), whereas the energy loss in 
forming the atom-pair defect at domain boundary is scaled as 0(1). From the analysis 
of the local energy of Fig. [3](b), the energy gain for each reconstructed 'red' atom is 
-E'gain ~ 0.6 eV (See the cases of the atom P or Q), while the energy loss of a 'green' 
atom of the defect is E\oss ~ 1-2 eV (See the cases of the atom R or S). The energy 
gain of helical domain is estimated to be 2i?gain per unit layer, since one unit layer of 
the wire contains two 'red' atoms. The energy loss at domain boundary, on the other 
hand, is estimated to be 2i?ioss, since the defect appears typically as an atom pair. The 
above estimation concludes that a helical domain should appear, if the wire length of 
the domain is enough long to satisfy Ldom > -^dom = (2-Eioss)/(2-E'gain) = 2 unit layers, 
which is consistent to the resultant domain structure of Fig. [2l 

We emphasis that the above scaling theory can be constructed, because the 
simulated wire is much longer than the critical length L^^^^^, unlike in the previous 
paper, in which the length is less than ten layers. [15] Although the atom movement 
from the inner shell into the wire surface is also found in a shorter wire of the 12-6-1 
structure, [15] the movement occurs within one atom line and the energy loss cannot be 
scaled as 0(1). 
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The present simulation gives two conclusive points for the formation of helical gold 
nanowire.; (i) Mechanism of atom insert on surface, which introduces helical domain 
and domain boundary with point defects, (ii) The simple energy scaling theory that 
explains the net energy gain. The above two points are universal and not dependent on 
the details of simulations. In general, the inserted atom can be supplied not only from 
the inner shell but also from the outer area, We speculate that another candidate 
for the source of the atom supply is the wire ends, if they are connected to electrode 
parts. 

5. Summary and future aspect 

In summary, domain boundary formation is explored for helical multishell gold nanowire. 
The simulations were carried out for nanowires longer than 10 nm by quantum 
mechanical molecular dynamics simulation. As results, the shear-like deformation on the 
(OOl)-type surface area introduces helical domains [15j and a domain boundary appears, 
typically with the supply of a defective atom pair from the inner shell, between two 
domains with the opposite shear directions. The mechanism is explained quantitatively 
by the analysis of local electronic structure and a general discussion is given by a simple 
energy scaling theory. 

As a future aspect, simulations should be carried out with larger samples that 
contain electrode parts, as pointed above. We note that the simulation with realistic 
electrodes is important also for the transport property, particularly among helical 
multishell nanowires and other nanowires thicker than the monoatomic chain, because 
of the possible interference effect at the connection part with electrodes. [l2l [6l [TJ HSj HI] 

A promising theoretical approach for larger quantum-mechanical simulations is 
'order-A^' method, in which the computational time is 'order-A^' or proportional to 
the system size A^. See articles cited in Ref. [33] or a recent journal volume that includes 
Ref. [16] and focuses on the order- A^ methods. 

References 

[1] Kondo Y and Takayanagi K 2000 Science 289 606 

[2] Oshima Y Kondo Y and Takayanagi K 2003 J. Electron Microsc. 52 49 
[3] Agrai't N Yeyati A L and van Ruitenbeek J M 2003 Phys. Rep. 377 81 
[4] Kondo Y and Takayanagi K 1997 Phys. Rev. Lett. 79 3455 
[5] Oshima Y Onga A and Takayanagi K 2003 Phys. Rev. Lett. 91 205503 
[6] Ono T and Hirose K 2005 Phys. Rev. Lett. 94 206806 

[7] Oshima Y Mouri K Hirayama H and Takayanagi K 2006 J. Phys. Soc. Jpn. 75 053705 
[8] Oshima Y Koizumi H Mouri K Hirayama H Takayanagi K and Kondo Y 2002 Phys. Rev. B 65 
121401(R) 

[9] Giilseren O Ercolessi F Tosatti E 1998 Phys. Rev. Lett. 80 3775 
[10] Tosatti E Prestipino S Kostlmeier S Dal Corso A and Di ToUa F D 2001 Science 291 288 
[11] Bilalbegovic G 2003 Vacuum 7 165 

[12] Senger R T Dag S and Ciraci S 2004 Phys. Rev. Lett. 93 196807 
[13] Lin J-S Ju S-P and Lee W-J 2005 Phys. Rev. B 72, 085448 



Domain boundary formation in helical multishell gold nanowire 



10 



Yang C-K 2004 Appl. Phys. Lett. 85 2923 

Iguchi Y Hoshi T and Fujiwara T 2007 Phys. Rev. Lett. 99 125507 

T. Fujiwara, T. Hoshi and S. Yamamoto, J. Phys.: Condens. Matter 20 294202 (2008). 

Amorim E P M and da Silva E Z 2008 Phys. Rev. Lett. 101 125502 

Iguchi Y 2007 D. Thesis (in Japanese) Univ. of Tokyo Tokyo 

Fedak D G and Gjostein N A 1967 Acta Metal. 15 827 

Binnig G Rohrer H Gerber Ch and Weibel E 1983 Surf. Sci. 131 L379 

Ho K-M and Bohnen K P 1987 Phys. Rev. Lett. 59 1833 

Van Hove M A Koestner R J Stair P C Biberian J P Kesmodel L L Bartos I Somorjai G A 1981 

Surf. Sci. 103 189 
Binnig G Rohrer Gerber C and StaU E 1984 Surf. Sci. 144 321 

Abernathy D L Mochrie S G J Zehner D M Griibel G and Gibbs D 1992 Phys. Rev. B 45 9272 
Jahns V Zehner D M Watson G M Gibbs D 1999 Surf. Sci. 430 55 
http: / /www. elses.jp/ 

Hoshi T and Fujiwara T 2000 J. Phys. Soc. Jpn. 69 3773 

Hoshi T and Fujiwara T 2003 J. Phys. Soc. Jpn. 72 2429 

Geshi M Hoshi T and Fujiwara T 2003 J. Phys. Soc. Jpn. 72 2880 

Takayama R Hoshi T and Fujiwara T 2004 J. Phys. Soc. Jpn. 73 1519 

Hoshi T Iguchi Y and Fujiwara T 2005 Phys. Rev. B72 075323 

Takayama R Hoshi T Sogabe T Zhang S-L and Fujiwara T 2006 Phys. Rev. B73 165108 
Hoshi T and Fujiwara T 2006 J. Phys: Condens. Matter. 18 10787 

Yamamoto S Hoshi T Sogabe T Zhang S-L and Fujiwara T 2008 J. Phys. Soc. Jpn. 77 114713 
Hoshi T and Fujiwara T 2009 J. Phys.: Condens. Matter. 21 064233 
Mehl M J and Papaconstantopoulos D A 1996 Phys. Rev. B 54 4519 

Kirchhoff F Mehl M J Papanicolaou N I Papaconstantopoulos D A and Khan F S 2001 Phys. Rev. 
B 63 195101 

Papaconstantopoulos D A and Mehl M J 2003 J. Phys: Condens. Matter. 15 R413 
Haftel M I and Gall K 2006 Phys. Rev. B 74 035420 

da Silva E Z da Silva A J R and Fazzio A 2001 Phys. Rev. Lett. 87 256102 

da Silva E Z Novaes F D da Silva A J R and Fazzio A 2004 Phys. Rev. B 69 115411 

Oshima Y Mouri K Hirayama H and Takayanagi K 2003 Surf. Sci. 531 209 

Kurui Y Oshima Y Okamoto M and Takayanagi K 2008 Phys. Rev. B 77 161403(R) 

Shinaoka H Hoshi T and Fujiwara T 2008 J. Phys. Soc. Jpn. 77 114712 



